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The Kibble-Zurek mechanism (KZM) predicts the density of topological defects produced in the 
dynamical processes of phase transitions in systems ranging from cosmology to condensed matter 
and quantum materials. The similarity between KZM and the Landau-Zener transition (LZT), 
which is a standard tool to describe the dynamics of some non-equilibrium physics in contemporary 
physics, is being extensively exploited. Here we demonstrate the equivalence between KZM in the 
Ising model and LZT in a superconducting qubit system. We develop a time-resolved approach to 
study quantum dynamics of LZT with nano-second resolution. By using this technique, we simulate 
the key features of KZM in the Ising model with LZT, e.g., the boundary between the adiabatic 
and impulse regions, the freeze-out phenomenon in the impulse region, especially, the scaling law 
of the excited state population as the square root of the quenching rate. Our results supply the 
experimental evidence of the close connection between KZM and LZT, two textbook paradigms to 
study the dynamics of the non-equilibrium phenomena. 

PACS numbers: 


Non-equilibrium phenomena at avoided level cross¬ 
ings play an essential role in many dynamical processes 
throughout physics and chemistry. A transition between 
energy levels at the avoided crossing is known as the 
Landau-Zener transition (LZT) [Il[2], which has served 
over decades as a textbook paradigm of quantum dynam¬ 
ics. LZT has recently been extensively studied [3] both 
theoretically and experimentally in, e.g., superconduct¬ 
ing qubits mz], spin-transistor [S] , and optical lattices 
[SHU]. On the other hand, quantum phase transition 
may also relate to avoided level crossings and it plays an 
important role in nature. Recently, an elegant theoreti¬ 
cal framework for understanding the dynamics of phase 
transition is provided by the Kibble-Zurek mechanism 
(KZM) [T3Mi)] . When the parameters of a quantum sys¬ 
tem that drive the quantum phase transition are varied 
in time causing the system to traverse the critical point, 
KZM predicts that the density of the defects produced 
in the processes follows a power law that scales with the 
square root of the speed at which the critical point is 
traversed. Due to its ubiquitous nature, this theory hnds 
applications in a wide variety of systems ranging from 
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cosmology to condensed matter and quantum materials 

[nMS]. 

The correspondence between LZT and KZM was first 
pointed out by Damski [211122] ■ It was shown that the 
dynamics of LZT can be intuitively described in terms 
of KZM of the topological defect production in nonequi¬ 
librium quantum phase transition. In order to model 
the dynamical process of LZT, a widely used picture is 
the adiabatic-impulse approximation (AIA), which was 
originally developed in KZM theory. The entire dynam¬ 
ical process can be divided into three regions: the adi¬ 
abatic, impulse, and adiabatic regions, as shown in Fig. 
la. The three regions are separated by two boundaries 
—vt and vt, where v is the quench rate and t is referred 
to as the freeze out time. Based on AIA, the dynam¬ 
ics of topological defect production in non-equilibrium 
phase transitions can be simulated with LZT. Recently 
this prediction was experimentally confirmed in an op¬ 
tical interferometer [23]. However, some key features in 
the correspondence between LZT and KZM, such as the 
freeze out time t and the adiabatic-impulse-adiabatic re¬ 
gions, have not been investigated experimentally. Most 
importantly, by studying the dynamical quantum phase 
transition in a quantum Ising chain, it is found that the 
average density of defects scales as the square root of 
the quenching rate [24| [2^. This universal scaling law of 
defect formulation as a function of quench time, which 
lies at the heart of KZM, lacks adequate experimental 
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FIG. 1: Energy level avoided crossing and experimen¬ 
tal procedure, (a) A typical energy structure (parame¬ 
terized by time) of a two-level system. The diabatic states 
|0), |1), and energy eigenstates lEip) are denoted in the plot, 
(b) A schematic of time profile of the experiment consisting 
of three parts. Dnring state preparation, a —7r/2j, (tt) pulse 
is applied to prepare the qubit in (|0) -I- |l))/\/2 (|1)). The 
Landau-Zener transition is realized by chirping the microwave 
frequency from Wi to ojf. The final state of the qubit is ob¬ 
tained by state tomography. 


evidence in LZT. 

In this paper, we use LZT in superconducting qubits 
to simulate KZM of the Ising model. We develop a time- 
resolved method to directly investigate the quantum dy¬ 
namics of LZT in the superconducting qubit. Using state 
tomography, we measure the time evolution of the popu¬ 
lation P_|_ of the instantaneous positive energy eigenstate 
[see Fig. la] for the entire LZT process. We find that 
P+ exhibits a rapid change near the center of the avoided 
crossing and varies gradually outside this region, reveal¬ 
ing the existence of the adiabatic and impulse regions. 
Moreover, the freeze-out behavior predicted by KZM has 
also been observed, and the boundary between the adia¬ 
batic and impulse region predicted by AIA is confirmed. 
We observe that the experimental simulated KZM of 
Ising model displays the theoretically predicted Kibble- 
Zurek scaling law. Therefore, we experimentally demon¬ 
strate the close connection between KZM and LZT, in 
particular, the presence of Kibble-Zurek scaling behavior 
in LZT. 

Results 

The equivalence between KZM and LZT. The Ising 
model is regarded as one of the two prototypical mod¬ 
els to understand quantum phase transitions | 2 I)] . After 
rescaling all the quantities to the dimensionless variables, 
we obtain the Ising model Hamiltonian 

N 

= + ( 1 ) 

n—1 

with the periodic boundary condition, where are the 
Pauli-matrices operators. Here N is the total number 
of spins and g is a dimensionless constant driving the 
phase transition. The ground state of Hj is a paramagnet 
for g ^ 1 and a ferromagnet for g <C I, and g = 1 
corresponds to the critical point. To study the dynamics 


of this model, we assume that the system evolves from 
time ti = —00 to = 0 , and takes a linear quench 
g{t < 0 ) = where 1 /tq provides a quench rate. 

By utilizing Jordan-Wigner transformation and 
Fourier transform, the Ising model can be simplified 
as a bunch of decoupled qubits with the Hamiltonian 

H'l = (lllk “ 5 ) j where the eigenenergy e/c = 

2^[g — cos kY + sin^ k and = - ,±(y — 5 )^ 

is the pseudomomentum. Here and hereafter we set 
the lattice constant a = 1. The density of defects re¬ 
sulting from the quantum quench has the expression 
N = '^kPk/N, where pk is the excitation probability 
corresponding to the pseudomomentum k. We consider 
a special case where the system undergoes slow evolu¬ 
tion with Tq ^ l/27r. Under this condition, it is safe 
to assume that only long wavelength modes are excited, 
i.e., k ^ 7 r/ 4 , and then pk « \uk(tf)\‘^, where {uk,Vk} 
are Bogoliubov modes governed by the following matrix 
equation 


= ( 2 ) 
dr \Uk J 2 \ 1 -TXk J \Uk J 

Here t = 4(t -|- Tq cos k) sin k is the normalized time, and 

Xfc = 1/4 tq sin^ k is the sweeping velocity. 

On the other hand, for a quantum two-level system 
with the diabatic basis , we consider the time-dependent 
Hamiltonian of a quantum two-level system in the dia¬ 
batic basis jo) and |1) 


H{t) 


1 ( e{t) A \ 

2 A -e(t) ) ’ 


( 3 ) 


where A is the tunneling amplitude and e(t) is the energy 
difference between the two diabatic basis. We mainly 
consider e(t) = vt with v being the speed of energy vari¬ 
ation. The time-dependent instantaneous energy eigen¬ 
states \E±{t)) are 


\E- (t)) = cos( 6 »/ 2 )| 0 ) -k sin( 6 l/ 2 )|I) . , 

\E+ (t)) = - sin(0/2)|O) + cos(d/2)|l) ’ 


where cos 9 = e(t)/H(t) with fl{t) = The 

instantaneous energy eigenvalues of H{t) are E±[t) = 
±|H(<), forming an avoided level crossing at t = 0 with 
a gap A. If the system initially in the ground state 
traverses the avoided crossing, the Landau-Zener theory 
gives the probability of the qubit occupying the excit¬ 
ing state as Plz ~ exp(— 7 rA^/ 2 Iir!). LZT has received 
tremendous attention since quantum two-level systems, 
i.e., qubits, are currently considered as the best building 
blocks of quantum information processors. 

With the substitution v/A^ — Xk and tA = r, it is 
found that the dynamics of the Ising model governed by 
Eq. (§ is the same as LZT physics contained in Eq. © 
up to a normalized tunneling amplitude A|25j. There¬ 
fore, we can use LZT to simulate KZM of the Ising model. 
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A key concept of KZM is AIA. Following the arguments 
in Ref. [211122] , we consider two nontrivial schemes to 
relate LZT with KZM. In scheme A, the system starts 
far away from the avoided crossing, corresponding to 
Ci —>■ —oo, and ends also far away from the avoided cross¬ 
ing, corresponding to e/ —>■ oo. The initial state |'0o) is 
the ground state of the Hamiltonian at time t = —oo. 
As shown in Fig. la, the evolution can be divided into 
three regions: there are two adiabatic regions t = (—oo, 
—t) and (i, oo), where almost no transition between 
the instantaneous energy eigenstates \E±{t)) occurs. On 
the contrary, in the impulse region [—i, t] transitions be¬ 
tween the states \E±{t)) could occur. Quantitatively, the 
boundaries separating the adiabatic and impuls e regions 
are determined by the freeze-out time t = ^/rQToJa, 
where a is a dimensionless parameter, tq = A/ v sets the 
scale of the quench time and tq = 1/A, respectively. The 
finite density of topological defects is caused by non- 
adiabatic evolution in the impulse region which 

equals the occupation probability of \E+). Therefore, the 
density of topological defects Dn in KZM corresponds to 
the transition probability P_|_ in LZT, i.e., = P.^.. 

In scheme B, the system starts from the center of the 
avoided crossing, i.e., Ci = 0, and evolves to the adia¬ 
batic region till far away from the avoided crossing. Sim¬ 
ilarly, there are two regions, an impulse region [0, i\ and 
an adiabatic {t, oo), can be defined. For both schemes, 
we directly measure the time-resolved of LZT in our 
experiment and then quantitatively compare the result 
with the prediction of KZM. 

The time-resolved LZT. We use superconducting 
qubits to investigate the dynamics of the Ising model. 
Two samples were studied: a superconducting phase 
qubit ( denoted as Qi with Ti = 113 ns, r 2 * = 93 ns) 
[6] and a 3D transmon ( denoted as Q 2 with Ti = 2.386 
/iS, T2 = 2.135 iis). Here Ti is the energy relaxation 
time from state |1) to state |0), and T 2 is the decoher¬ 
ence time including contributions from both relaxation 
and dephasing. Because neither the phase qubit nor the 
transmon qubit possess an intrinsic avoided level cross¬ 
ing, we use a coherent microwave field to generate an ad¬ 
justable effective avoided energy level crossing |27|. The 
position and the tunneling amplitude (A) of the avoided 
crossing are determined by the frequency and amplitude 
of the microwave field, respectively. With this flexibil¬ 
ity and controllability, instead of sweeping the flux bias, 
we chirp the microwave frequency ut while keeping the 
qubit frequency wio constant to realize LZT (see Meth¬ 
ods) . The profiles of the control and measurement pulses 
are illustrated in Fig. lb. 

In our experiment, scheme A is a good approximation 
to LZT from —00 to -l-oo because ei/27r = —200 MHz and 
A/27r = 20 MHz, resulting je^j/A = 10. The initial state 
is prepared in |1) « \E-) at ei/27r = —200 MHz. This is 
followed immediately by chirping e to e/ with a constant 
speed V = [cf — where t^z is the duration of the 

chirping operation. At the end of e chirping, we perform 
state tomography to determine the density matrix of the 
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FIG. 2: The values of {ox.z) and P+ as a function of 

e//27r and Ilz- Here = —200 MHz and A/27r = 20 

MHz. The range of e//27r is from —200 MHz to 400 MHz. The 
LZ duration Ilz is from 1 ns to 120 ns. (a)-(c) ((d)-(f)) are 
the experimental (numerically simulated) results. 


qubit p = + a ■ (a )) by measuring the expectation 

values {<Jx,y,z), where are the Pauli matrices. We 

then varied tzz from 1 ns to 120 ns, and e//27r from 
—200 MHz to -1-400 MHz to obtain a complete set of 
experimental data. Since the avoided crossing centers at 
e = 0, our experiments effectively cover the dynamical 
evolution of the system from e<C—Atoef^A. 

By converting the density matrix to the time- 
dependent basis \E±), we obtain the 

P+{t) = {E+{t)\p\E+{t)) 

1 (5) 

= 2 ( 1 “ (cr^) cos d- (cr:j,)sm6l), 

which shows that only {az) and (ct^) contribute to P+- 
The measured (cr^), (a^) and P+ for qubit Q 2 are plotted 
in Fig. 2a, b and c, respectively. Shown in Fig. 2d-f are 
the results of numerical simulation obtained by solving 
the master equations (see Methods), where important 
system parameters, such as the relaxation and the de¬ 
phasing times are determined from the pump-decay and 
the Ramsey fringe measurements. The good agreement 
between the experimental and the simulated results indi¬ 
cates that all essential aspects of our experiment are well 
controlled and understood. 

The adiabatic and impulse regions. The time re¬ 
solved quantum dynamics of P+ (t) described above can 
be investigated by measuring {<Jx,z) with nano-second 
time resolution. Shown in Fig. 3a are examples of P+ 
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FIG. 3: Population P+ as a function of the normalized 
time t/t and the comparison of T>n and (a) the 

evolution starting from t = —oo with ei/2'n = —200 MHz 
and e//27r = 200 MHz. (b) the evolution starting from t = 0 
with ei/27r = 0 and e//27r = 400 MHz. Different LZ durations 
tnz = 10 ns (red circle), 20 ns (magenta square), 40 ns (blue 
triangle), 80 ns (green diamond), are used to produce differ¬ 
ent LZT speed. The symbols (solid lines) are experimental 
(numerical) results. The red translucent (clear) regions mark 
the impulse (adiabatic) regions, while the boundary locates 
on ±t. The error bars are smaller than the sizes of the sym¬ 
bols. (c) The comparison of topological defects density 
in KZM theory and P+ in LZT with ei/27r = 0. The blue 
symbols (green solid lines) are the experimental (numerical) 
results. The red dashed line shows the density predicted 
in KZM with a = 0.784 as the best fit. 


for qubit Q 2 as a function of evolution time for various 
LZT duration time t^z- In order to compare with AIA, 
we normalize the evolution time by the freeze-out time 
i = ^Xjav. The value of ot used here is 7r/2, which is the 
same as that of AIA in this scheme Enna. In all cases 
P+ changes rapidly near the center of the avoided cross¬ 
ing and varies slowly outside the central region. This 
is the clear experimental evidence supporting the physi¬ 
cal picture of AIA. Moreover, the boundaries between the 
adiabatic and impulse regions are demarcated by ±i with 
no fitting parameters, confirming the validity of AIA. 

In scheme B, we investigate LZT by starting from the 
center of the avoided crossing (i.e., Ci = 0). The system 
is initialized in the lower energy eigenstate at = 0 
with a proper resonant 7r/2 pulse [251130) . Then a time 
sequence similar to that in scheme A is applied. Here, 
e//27r ranges from 0 to 400 MHz. The gap size is fixed at 
A/27r = 20 MHz, resulting a maximal e//A = 20. Shown 
in Fig. 3b are examples of measured P+ as a function of 
the evolution time. In this case, a = T^jA according to 
Ref. [211 Eg. Similar adiabatic and impulse regions are 
observed with the boundary at t = t, strongly supporting 
AIA. 

The freeze-out phenomenon. Another interesting 
problem is whether one can observe directly the predicted 
state freeze-out phenomenon in the impulse region [—t, i]. 
According to KZM, although P± of the time-dependent 
basis states \E±{t)) change rapidly in the impulse re¬ 
gion, the probability amplitudes (thus ("d^)) in the time- 
independent basis {|0), 11)} should be frozen out. To see 
whether this is indeed the case, we plot the measured 
(o'z.x.y) of the qubit Q 2 in Fig. 4a-c (see Methods). The 
line represents the freeze-out time t = Xj^fon) is also 



FIG. 4: State freeze-out phenomena. (a)-(c) ((d)- 
(f)) are the experimental observation (numerical simulation) 
of the state freeze-out phenomena of the expectation values 
{(yz,x,y) in LZT with ei/27r = 0. The white solid line marks 
the freeze-out time t in KZM with a = 7r/4. 


shown in the plot. Here we use the theoretical KZM 
value a = 7r/4 because the total LZT duration is shorter 
than 40 ns and the effect of decoherence may not have sig¬ 
nificant effects on the result. It can be seen that {<Jz^x,y) 
change slowly in the impulse region, indicating that the 
state of the qubit is nearly frozen. In order to compare 
with the experimental data, we present the numerically 
simulated results by solving master equations without 
adjustable parameters for the evolution of LZT in Fig. 
4d-f. The good agreement between the simulation and 
the experimental results supports the observation of the 
state freeze-out phenomenon and confirms the validity of 
AIA. 

One of the key features of the correspondence between 
KZM and LZT is ~ P+. In Fig. 3c we plot as a 
function of t^z for e//27r = 200 MHz thus |e/|/A = 10. 
In order to compare with the theory, t^z is expressed in 
terms of tqItq. It is found that P+ follows quite well 
with the behavior of the topological defects density 
predicted in Ref. [211122]. Here, a = 0.784 is obtained 
from the best fit which is within 0.2% of KZM predicted 
value 7r/4. The excellent agreement between the exper¬ 
imental results of LZT and the theory of KZM provides 
strong support to the conjecture that the dynamics of 
the Landau-Zener model can be accurately described in 
terms of the Kibble-Zurek theory of the topological de¬ 
fect production in nonequilibrium phase transitions and 
vice versa. 

The scaling law. We now address the simulation of 
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the scaling law predicted by KZM for the Ising model. 
By choosing small quenching rates 1 /tq, all the rele¬ 
vant physics described by Eq. (2) happens in the long 
wavelength limit fc <C tt. In experiments, we choose a 
cutoff /cc/tt and thus Uc/A^ = Xkc to ensure that LZT 
probability can be neglected for |A:| > k^- For each Tq, 
we choose Ni~ different quasimomentum k equally dis¬ 
tributed in [—kc,kc], and measure the corresponding ex¬ 
citation probability P+{k,TQ). Then the average density 
of defects can be expressed as 




irNk 



( 6 ) 


where the last equation is given by KZM theory. Stimu¬ 
lated by this prediction, we plot experimentally measured 
A/’(tq) vs. 1/y^ qubits Qi (red squares) and 

Q 2 (blue squares) in Fig. 5, where Nk = 127 for each 
Tq, kc/iT = 0.2 and 1 /tq < 0.01. A striking feature 
is that M{tq) shows a very good linear relation with 


1 


Q- 


By fitting the line to a general linear function 


— Nq-\- 15!^Tq^ we obtain the offset Ag and slope 
P as summarized in Table 1. It is interesting that with the 
increasing of the decoherence time, the slope increases 
while the offset decreases. In order to confirm our obser¬ 
vation, we did numerical simulations using Ti and T 2 of 
Qii Q 2 , and infinite, shown in Fig. 5. Since there is no 
adjustable parameter, the agreement between the experi¬ 
mental data and numerical simulation results are remark¬ 
able. When the decoherence time goes to infinite, the off¬ 
set tends to zero and the slope P ~ 0.106 ± 0 . 002 , which 
is very close to the theory predicted value [25] . 

Therefore, it is evident that LZT in qubit exhibits same 
scaling behavior of KZM of the Ising model although de¬ 
coherence will quantitatively modify the value of some 
parameters. 

Discussion 

In our experiment, the hallmark features of the Ising 
model predicted by KZM, such as the existence of the adi¬ 
abatic and impulse regions, the freeze-out phenomenon, 
and scaling law, were observed. The experimental obser¬ 
vations of the first two are in good agreement with the 
theoretical results because in the experiment the time of 
evolution is much less than Ti and T| of the qubit Q 2 - 

We here make several comments on our results of the 
scaling law. (i) In the absence of decoherence, our nu¬ 
merical results show that the slope P ~ 0.106, while the 
theory predicts P ~ 0.1125. The discrepancy between 
the numerical and the theoretical results is due to the 
fact that in the numerical simulation, both and e/ 
are finite, the LZ transition happens in a finite range, 
i.e., e//A ~ 10. In theory, however, P is assuming 
e//A —>■ 00 . To confirm it, we increased e//A in nu¬ 
merical simulation, and indeed find that P approaches 
\/2'/2tt asymptotically, (ii) The offset Nq in the defect 
density N is observed for both qubits in Fig. 5. In addi¬ 
tion, the stronger the decoherence, the larger the offset. 



FIG._5: 


Q- 


The scaling behavior of A as a function of 

The red (blue) squares represent the experimental 


data measured in Qi (Q 2 ). The magenta, blue solid lines 
are the numerical simulation of the master equation with the 
decoherence of the phase qnbit and 3D transmon, respectively. 
The green solid lines are the simnlated results with infinite Ti 
and T 2 ■ 


It was found in previous studies |3T1132| that decoherence 
tends to increase the density of defects. As the decoher¬ 
ence rate increases, more defects are generated in the 
process of LZ transition. This is confirmed in our exper¬ 
iments, where it is found that the defect density N as 
well as the offset Ag for qubit Qi is greater than that 
for qubit Q 2 - If the coherence time increase further, the 
offset will gradually approach to zero, as verified by the 
simulation 3 in Fig. 5. (iii) It is known that the theoret¬ 
ical prediction of Ref. |25| is obtained in the continuum 
limit. So a further question is whether the finite size ef¬ 
fect is important in our experiments. If we consider a 
quantum Ising spin chain with a finite length A, there 
will be energy splitting in the energy spectrum. In this 
case, if the chirping speed is too slow compared with the 
energy splitting caused by the finite size, then the predic¬ 
tion of Ref. [25] would not be valid. In other words, the 
finite size effect sets a lower bound on the chirping speed. 
To be more precise, let us consider the energy splitting 
near zero energy where k = ± 5 ^. Based on Eq. (j^, we 

require Xk=±n/N > 1, corresponding to > ^. In 

our experiments, we choose A = Nkir/kc — 1000, which 
requires y^I/tq > 0.006, a condition that is very well 
satisfied in our experiments. Thus the effect of finite size 
is negligible. 

In conclusion, using linear chirps of microwave field 
and exploring the correspondence between the KZM of 
topological defects production and LZT, we simulated 
the KZM of the Ising model with a superconducting qubit 
system. All important predictions of KZM for the Ising 
model, such as the existence of adiabatic and impulse 
regions, the freeze-out phenomenon, and especially the 
scaling law have been clearly demonstrated. The ob- 
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Samples 

Ti 

n 


No 

Qi 


113 ns 

93 ns 

0.068 ± 0.002 

0.0091 ±0.0001 

Q2 


2.386 ps 

2.135 ps 

0.088 ±0.002 

0.0048 ± 0.0001 

Simu. 

1 

113 ns 

93 ns 

0.070 ±0.001 

0.0087 ±0.0001 

Simu. 

2 

2.386 ps, 

2.135 ps 

0.090 ±0.002 

0.0048 ± 0.0001 

Simu. 

3 

00 

00 

0.106 ±0.002 

-0.0001 ±0.0001 


TABLE I: The offset No and slope /3 extracted from the 
experimental (Qi, Q 2 ) and simulated results (Simu. 1, Simu. 
2, Simu. 3). 


served scaling behavior in the presence of decoherence 
sheds new light on the investigation of the effects of de¬ 
coherence on KZM of non-equilibrium quantum phase 
transitions. 

Methods 

Chirp-LZT operation. In order to perform the LZT, 
we chirp the microwave frequency |29[ I30| instead of 
sweeping the flux bias of the qubit. Concretely, with the 
qubit dc-biased at a fixed flux $, we chirp the microwave 
frequency from uji to w/, corresponding to the change of e 
from Ci to ef. Note that in our experiment, the energy dif¬ 
ference is Cij = h{ujQi — ujij), with a chirping microwave 
frequency uJij = Wgi + f Here we assume h = 1. 
If we set the original microwave frequency as wqi, then 
the chirped frequency is ^ ■ Therefore, the relationship 
between the chirped microwave frequency (5^^- ^ and the 
diabatic energy difference Cij is given by Cij = —5^. ^. 

To chirp the microwave frequency, we apply modula¬ 
tion signals from a Tektronix AWG70002 to the IF (in¬ 
termediate frequency) ports of a IQ mixer. Considering 
the original microwave waveform as A^sinwot, the mod¬ 
ulation signals applied to the I and Q ports of the IQ 
mixer as cos Suit and sin Sujt, respectively. In this way, the 
modulated microwave at the output port of the mixer is 
Aj, sinwg^ cos -I- coswotsind^^t = sin(cLio + Sui)t, 
and the microwave frequency is chirped by tuning Sui- 
From calibration, we find that the power of luq tone is 
at least 50 dB lower than that of ujo + Sui, indicating the 
negligible effect of loq on the qubit in the chirp operation. 

The Chirp-LZT method provides us several advantages 
in performing LZT. First of all, all parameters of the 
qubit, such as Ti, T 2 and the coupling strength between 
the qubit and the external driven filed, are fixed dur¬ 
ing the measurements. Second, the end points of dia- 
batical energy sweep and e/ would not be limited by 


the avoided level crossings resulting from the coupling 
between the qubit and microscopic two level systems 
(TLSs) usually presented in superconducting qubits be¬ 
cause the microwave couples much weakly to TLSs than 
to the qubit. Third, it is easy to control the chirping 
velocity and the tunnel splitting A by controlling the 
frequency and power of the microwave. In addition, such 
method provides a useful tool for systems without natural 
avoided crossings in their energy diagram, such as trans- 
mons, to perform LZT and other similar experiments. 
Solution of the master equation. The numerical 
results are obtained by solving master equations. The 
quantum dynamics of the superconducting qubits is de¬ 
scribed by the master equations of the time evolution 
of the density matrix p including the effects of dissipa¬ 
tion: p = — r[p], where H is the Hamiltonian of 

the system described by Eq. (3). The second term of 
this equation, r[p], describes the effects of decoherence 
on the evolution phenomenologically. Setting h = 1, the 
master equation can be rewritten as 


{ Poo = |^(pio — Poi) + r ipii 

Pii = |A(poi - Pio) - Tipii 

Poi = |[A(pii - Poo) + 2epoi] - 7Poi ^ 
Pw = iiA(poo - Pii) - 2epio] - 7Pio- 

Here Fi = 1/Ti is the energy relaxation rate, and 7 = 
I/T 2 is the decoherence rate.The relationship between 
the density matrix and the qubit state expectation values 
are given by 


( Poo — (1 + (o'z))/2 

1 Pii = (1 - (o-2)V2 

I ^e(poi) = (cra;V2 
I Im{poi) = -{uy)/2 


( 8 ) 


The numerical simulations in Fig. 2-4 are straightfor¬ 
wardly obtained by solving Eq. 0 . The simulations 
(solid lines) in Fig. 5 are obtained by mapping P+ as 
defined in Eq. ([^ to the defect density M. By using 
ri/A^ = Xk for each chirping velocity v, we can find the 
corresponding momentum k with Tq fixed. In this way, 
we obtain P+ik,rl^) for different momentum k. Then 
A/’(tq) is computed from P+ using Eq. (j^. 
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